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ABSTRACT 

This work investigates the main mechanism(s) that regulate the specific star formation 
rate (SSFR) in nearby galaxies, cross-correlating two proxies of this quantity - the 
equivalent width of the Ha line and the {u — r) colour - with other physical properties 
(mass, metallicity, environment, morphology, and the presence of close companions) 
in a sample of 82500 galaxies extracted from the Sloan Digital Sky Survey (SDSS). 
The existence of a relatively tight ‘ageing sequence’ in the colour-equivalent width 
plane favours a scenario where the secular conversion of gas into stars (i.e. ^nature') is 
the main physical driver of the instantaneous SSFR and the gradual transition from 
a ‘chemically primitive’ (metal-poor and intensely star-forming) state to a ‘chemi¬ 
cally evolved’ (metal-rich and passively evolving) system. Nevertheless, environmental 
factors (i.e. ^nurture’’) are also important. In the field, galaxies may be temporarily 
affected by discrete ‘quenching’ and ‘rejuvenation’ episodes, but such events show lit¬ 
tle statistical significance in a probabilistic sense, and we find no evidence that galaxy 
interactions are, on average, a dominant driver of star formation. Although visually 
classified mergers tend to display systematically higher EW(Ha) and bluer {u — r) 
colours for a given luminosity, most galaxies with high SSFR have uncertain mor¬ 
phologies, which could be due to either internal or external processes. Field galaxies 
of early and late morphological types are consistent with the gradual ‘ageing’ scenario, 
with no obvious signatures of a sudden decrease in their SSFR. In contrast, star for¬ 
mation is significantly reduced and sometimes completely quenched on a short time 
scale in dense environments, where many objects are found on a ‘quenched sequence’ 
in the colour-equivalent width plane. 

Key words: galaxies: star formation - galaxies: emission lines - galaxies: interactions 
- galaxies: abundances 


1 INTRODUCTION 

Galaxies are known to evolve through a mixture of intrinsic 
(‘nature’) processes, such as gas accretion, conversion into 
stars, and chemical enrichment, as well as through interac¬ 
tions (‘nurture’) with other galaxies. The observed distribu¬ 
tion of their physical properties (e.g. stellar and gas mass, 
metallicity, SSFR, luminosities and colours) is set by a more 
or less complex combination of internal and external pro¬ 
cesses, but their relative importance is still a matter of de¬ 
bate. In this work, we are interested in the mechanisms that 
regulate the (specific) star formation rate in nearby galax- 
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ies, and, more specifically, the role of dense environments 
and close interactions on quenching and/or triggering star 
formation. 

Over the last decades, many studies have found that 
star-forming galaxies (the so-called ‘blue cloud’) are neatly 
segregated from passively-evolving systems (the ‘red se¬ 
quence’) in a colour-magnitude diagram (e.g. Tully et al. 
1982; Strateva et al. 2001; Baldry et al. 2004, 2006). Objects 
in the intermediate region, known as the ‘green valley’, are 
typically considered a transition population (e.g. Bell et al. 
2004; Faber et al. 2007; Martin et al. 2007; Schiminovich 
et al. 2007; Wyder et al. 2007; Mendez et al. 2011; Gongalves 
et al. 2012). Many of them display signatures of active galac¬ 
tic nuclei (AGN) in the optical spectrum (e.g. Salim et al. 
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2007; Schawinski et al. 2007, 2010), and their colours are 
often interpreted as evidence for a putative ‘quenching’ pro¬ 
cess that would quickly interrupt any star formation activity. 

In this work, we will use the word ‘quenching’ in the 
sense of a discrete event that instantaneously (or in a very 
short time scale) truncates the star formation activity in a 
given galaxy. In contrast, we will use the term ‘ageing’ to 
denote a gradual increase in the mass-weighted average age 
of the stellar population 


r* (t) = t — 
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without any sudden change in its star formation rate (SFR) 
'tp(t). With this definition, it can be easily shown that galax¬ 
ies become ‘older’ as they evolve unless their stellar mass 
increases faster than exponentially (e.g. during a strong 
starburst). Therefore, all galaxies are subject to ‘ageing’, 
whereas they may or may not undergo ‘quenching’ or ‘reju¬ 
venation burst’ episodes. 

It has long been known that galaxy colours are strongly 
correlated with environment, suggesting that ‘quenching’ 
might be a ‘nurture’ process. The location of the red se¬ 
quence is fairly independent on galaxy density (Sandage 
& Visvanathan 1978), and the position of the blue cloud 
in the colour-magnitude diagram varies only weakly (e.g. 
Balogh et al. 2004). Most galaxies in the blue cloud display 
similar SSFRs, varying relatively slowly with mass (Noeske 
et al. 2007; Speagle et al. 2014), and it has been shown (e.g. 
Balogh et al. 2004; Tanaka et al. 2004; Park et al. 2007; 
Wijesinghe et al. 2012) that the distribution of Ha equiva¬ 
lent widths for the star-forming population, selected on line 
intensity and/or colour, is not a strong function of local 
density. However, the fraction of red and blue galaxies is a 
strong function of the environment. 

One possible interpretation would be that all galax¬ 
ies evolve along a ‘star-forming sequence’ until a quenching 
event would rapidly drive them through the green valley to¬ 
wards the red sequence (see e.g. Peng et al. 2010). According 
to these authors, there are two different quenching mecha¬ 
nisms, one related to mass (i.e. ‘nature’) and one related 
to environment (‘nurture’). As galaxies grow in mass, some 
unidentified process would shut down star formation with a 
probability that scales linearly with the instantaneous star 
formation rate (as opposed to e.g. a fixed mass threshold). 
‘Environment quenching’ would be driven by galactic over¬ 
density, and it is found to be roughly independent of stellar 
mass or cosmic epoch up to 2 ; 1, consistent with the infall 

of the galaxy into a larger dark matter halo. Both mecha¬ 
nisms are independent, and their total effect on the fraction 
of red galaxies is fully separable. 

This kind of scenarios implicitly assume that star¬ 
forming and passive galaxies form two distinct groups. How¬ 
ever, some studies question the very existence of such a bi¬ 
modality in the galaxy population. Applying data mining 
techniques, Ascasibar & Sanchez Almeida (2011) found that 
almost all SDSS galaxies are distributed along a well-defined 
curve in the multidimensional space defined by their spectra. 
Only optically-bright active galaxies appear as an indepen¬ 
dent, roughly orthogonal branch that intersects the ‘main 
sequence’ exactly at the point of the transition between star¬ 
forming and passive systems. The location of such transition 


is thus well defined, but there is no apparent gap, as far as 
the optical spectra are concerned. 

As recently pointed out by Schawinski et al. (2014), 
galaxy morphology may be an important piece of the puzzle. 
While early-type galaxies display signatures of rapid quench¬ 
ing, leading to a fast transition through the green valley, 
late-type galaxies are consistent with slowly-declining star 
formation rates, with typical time scales of the order of sev¬ 
eral Gyr. They do not separate into a blue cloud and a red 
sequence, but rather span almost the entire (u — r) colour 
range without any gap or valley in between. The apparent 
bimodality in the colour-magnitude diagram arises from the 
superposition of both morphological populations. 

Schawinski et al. (2014) identify early-type galaxies 
with merger-induced (i.e. nurture) quenching, whereas they 
advocate for a discrete quenching event in the distant past 
of late-type galaxies, perhaps associated with reaching a 
certain halo mass (i.e. nature), that would have shut off 
the infall of new cold gas. However, the instantaneous SFR 
would not drop immediately, but gradually decrease until 
the cold gas reservoir is exhausted. Therefore, this process 
would qualify as ‘ageing’ rather than ‘quenching’ according 
to our terminology. 

Here we would like to statistically quantify the secu¬ 
lar evolution of the SSFR as well as the relevance of envi¬ 
ronment and interactions. Galaxies in dense environments 
are subject to a plethora of processes, such as tidal forces, 
strangulation of the gas supply, ram-pressure stripping, or 
galaxy harassment (see e.g. Boselli & Gavazzi 2006) that 
would suppress the star formation rate, although the exact 
time scale is still debated (see e.g. Wijesinghe et al. 2012, 
and references therein). 

On the other hand, close interactions with neighbour¬ 
ing galaxies can lead to gas compression and trigger star 
formation (see e.g. Li et al. 2008; Ellison et al. 2008; Ideue 
et al. 2012). Extreme objects such as ultraluminous infrared 
galaxies are always found to be strongly interacting systems 
(Sanders & Mirabel 1996; Surace et al. 1998). However, the 
mechanisms that set the star formation rate at the low-mass 
end are yet unclear. Some of these systems show very high 
SSFRs, often interpreted as massive ‘rejuvenation’ bursts 
of star formation (e.g. Heckman et al. 1998), but there is 
little evidence that these episodes are triggered by interac¬ 
tions with nearby neighbours (see e.g. Fujita 1998; Martig 
& Bournaud 2008; Lamastra et al. 2013; Lanz et al. 2013, 
and references therein). HII galaxies represent the most ex¬ 
treme objects in this mass range. They are, in general, blue 
compact dwarf (BGD) galaxies showing very strong and nar¬ 
row Balmer emission lines with large equivalent widths (Sar¬ 
gent & Searle 1970), and they have the lowest metal content 
of any star-forming galaxy (Searle & Sargent 1972; Rosa- 
Gonzalez et al. 2007). Observations show that the most lu¬ 
minous HII galaxies tend to live in very low density en¬ 
vironments without obvious companions (see e.g. Telles & 
Terlevich 1995; Vilchez 1995), although they show irregular 
morphologies and large velocity dispersions, consistent with 
a ‘nurture’ origin (Telles et al. 1997). 

The main goal of the present work is to investigate the 
physical properties of star-forming galaxies in the local uni¬ 
verse in order to disentangle the contribution of nature and 
nurture to galaxy ageing, quenching, and rejuvenation. In 
particular, we aim to explore the role of stellar mass, local 


© 2015 RAS, MNRAS 000, 1-18 



SSFR in SDSS galaxies: nature or nurture? 3 


environment, close interactions, and metallicity in setting 
the SSFR. Although there is a vast amount of literature on 
this topic, this work presents one of the most comprehen¬ 
sive studies of all the relevant galaxy properties to date. 
It has also sufficient statistics to provide a fair represen¬ 
tation of very faint objects, extreme values of the SSFR, 
and evolutionary phases that last for relatively short times. 
We have chosen to use purely observable quantities rather 
than model-dependent inferences, and we propose new trac¬ 
ers of galactic overdensity and close interactions. Most im¬ 
portantly, we propose a physical interpretation of our results 
that is somewhat different from the currently most widely 
accepted scenarios. 

Our selection of a suitable galaxy sample and its phys¬ 
ical characterization are described in Section 2. The distri¬ 
bution of the different observables we have considered, their 
correlations, and the implications regarding the mechanism 
that regulates star formation are discussed in Section 3. A 
theoretical interpretation of our results is discussed in Sec¬ 
tion 4. Finally, our main conclusions are summarized in Sec¬ 
tion 5. 

2 OBSERVATIONAL DATA 
2.1 Sample definition 

In order to proceed with the analysis we require a large 
sample of galaxies (hereafter main) that will be used to de¬ 
fine neighbours and environment, as well as to select the set 
of subsamples whose physical properties will be discussed 
throughout the paper. 

All data used in the present analysis are obtained 
from the SDSS Data Release 7 (Abazajian et al. 2009) 
database. Only objects within the completeness threshold 
of nir < 17.77 and catalogued as galaxies by the spectro¬ 
scopic pipeline (entries listed in table SpecObj^) were cho¬ 
sen to become part of the main sample. Line measurements 
are obtained from table SpecLine, and absolute magnitudes 
in the r band, Mr, have been computed from the apparent 
model magnitude^ listed in the photometric catalogue (Pho- 
toObjAll table) using the spectroscopic redshift distance. 

In addition, we impose the following criteria to select 
our target galaxies from the main sample: 

(i) 132.0 < RA < 230.0 and 1.0 < DEC < 56.0 

(ii) TTlr < 17.5 

(iii) 0.02 < 2 : < 0.07 

As can be readily seen in Figure 1, the region of the sky 
from which the star-forming subsample is selected (white 
square) leaves at least 2° (~ 3 and 10.5 Mpc for z = 0.02 
and 0.07, respectively) with respect to the main sample foot¬ 
print, in order to minimize boundary effects when studying 
interactions and environment. The restriction in apparent 
magnitude (m-r < 17.5) ensures that every object with a 
similar-mass companion would be properly identified. The 
upper redshift considered roughly corresponds to an abso¬ 
lute magnitude Mr < —20.0, which is approximately the 
threshold below which most galaxies are actively forming 

^ http://cas.sdss.org/dr7/en/help/browser/browser.asp 
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Figure 1. Spatial distribution of the main sample of galaxies, 
and region of the sky from which the star-forming subsample is 
selected. 
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Figure 2. Distribution of galaxies with S/N > 2 measurements 
of the emission lines in the BPT diagram. A dotted line shows 
the Kewley et al. (2006) criterion used to separate star-forming 
galaxies from AGN. 


stars, and the lower redshift cut prevents selection effects 
due to local structure. 

We then classify our target galaxies into four different 
subsamples: 

(i) Galaxies with signal-to-noise ratio S/N > 2 in the 
Ha, H/3, [OIII], and [Nil] emission lines, classified as AGN 
according to the Kewley et al. (2006) criterion based on 
their position on the BPT diagnostic diagram (Baldwin et al. 
1981, see Figure 2) are labelled as agn and will not be con¬ 
sidered further. 

(ii) Galaxies with S/N > 2 in all the above lines, classified 
as star-forming according to Kewley et al. (2006), will be 
referred as such hereafter. For these systems, the imposed 
S/N ensures that the quality of the spectra is sufficient to 
estimate the gas-phase metallicity from the 03N2 ratio. 

(iii) Galaxies where Ha is observed in absorption with 
S/N > 2 will be referred to as passive. 

(iv) All other objects are classified as intermediate. 
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2.2 Physical characterization 


In order to characterize our galaxy samples in a model- 
independent fashion we will use direct observables as proxies 
for their physical properties. 

First and foremost, we will use two different indicators 
of the SSFR, which is the main focus of the present study. On 
the one hand, we will consider the equivalent width of the Ha 
emission line, EW(Ha), where the line intensity traces the 
presence of young stars while the underlying continuum is 
set by the whole star formation history of the galaxy (see e.g. 
Terlevich et al. 2004). On the other hand, we will also use the 
(u —r) colour obtained from the model magnitudes, sensitive 
to the presence of massive, short-lived stars that dominate 
the bluer spectral bands. The main differences between both 
observables are that EW(Ha) is subject to aperture effects 
(see below) and that it traces a younger stellar population 
than (u — r). Therefore, these two indicators probe the SSFR 
on different spatial and temporal scales. 

In order to describe the intrinsic properties (i.e. ‘na¬ 
ture’) of our objects, throughout the present work we will 
use the absolute magnitude in the r band, Mr, as a proxy 
for stellar mass, and we will use both terms interchange¬ 
ably. To account for their evolutionary state, from the point 
of view of chemical enrichment, we will use the observational 
indicator 


03N2 = log 
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that is well known to decrease monotonically with the gas- 
phase metallicity (Alloin et al. 1979) and thus with stellar- 
to-gas fraction (Searle & Sargent 1972; Edmunds 1990; As- 
casibar et al. 2014). 

Concerning the effects of ‘nurture’ on the SSFR, we 
will quantify the environment in terms of the local galaxy 
overdensity, estimated from the projected distance to the 
fifth nearest neighbour in the main sample. We propose the 
relative distance between the first and second neighbours as 
a simple indicator for the presence of close companions, and 
we will interpret galaxy morphology as a possible tracer of 
galaxy-galaxy interactions. 

The definition of distance is a key aspect in the char¬ 
acterization of close interactions and environment. Spectro¬ 
scopic redshifts are available for every object in our sample 
which permits distance calculations in 3D space. However, 
redshift measurements are strongly affected by the peculiar 
velocities of galaxies (especially important in denser environ¬ 
ments such as galaxy clusters). As our galaxy sample covers 
all kind of environments, this effect will act irregularly and 
introduce a bias. In order to avoid such problems, we base 
our clustering calculations on projected distance measure¬ 
ments, as it is widely done in the literature (see e.g. Kaiser 
1987; Cooper et al. 2005; Muldrew et al. 2012, and references 
therein). 

For every target galaxy with redshift we select po¬ 
tential neighbours from the main sample within a range 
2gai ± 0.005 (which corresponds to 15 Mpc or 1500 
km/s velocity difference). Then, the distance to each object 
in this redshift shell is calculated as the great-circle distance, 
i.e. projected onto the surface of a sphere with radius dA{z), 
where 


dA{z) 


dL{z) 
{l + zY 


(3) 


is the angular diameter distance, and dL{z) corresponds to 
the luminosity distance at redshift 2 . 

Since our main sample is magnitude-limited, the num¬ 
ber density of galaxies decreases with redshift as the faintest 
objects drop out from the sample. Hence, the average dis¬ 
tances to the neighbours increase with z. We circumvent 
this problem by normalizing such distances to a character¬ 
istic radius 


r{z) 


dAjz) 
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that is proportional to the average distance between galax¬ 
ies, assuming they were uniformly distributed over the solid 
angle AO within the redshift shell 2 ± 0.005 (i.e. Trr^ oc 
AOd^/Agai). A^gai refers to the number of galaxies within 
the redshift shell under consideration. Then, we will use 
the normalized distance to the fifth neighbour, Rs/f, as a 
redshift-independent estimator of the local galaxy density 
(i.e. an indicator of the environment where the galaxy lives 
in). Galaxies in voids display values of the order of a few, 
whereas cluster galaxies feature values well below unity. 

In order to test for the presence of close companions, 
we chose to use the ratio between the squared projected 
distances to the first and second neighbours. For a uniform 
random distribution in two dimensions, R 1 /R 2 should be 
uniformly distributed between 0 and 1, irrespective of the 
local galaxy density. A low value of Ri is not in itself in¬ 
dicative of interactions, and it is trivially obtained e.g. for 
all cluster galaxies. In contrast, very low values of R 1 /R 2 
may mean (in a probabilistic sense) that the first neighbour 
is much closer than expected for a purely random distribu¬ 
tion, and thus they trace the fraction of close galaxy pairs 
in a statistical sample. 

Finally, the morphology of a galaxy is also a powerful 
indicator of its dynamical and merger history. Thanks to the 
Galaxy Zoo 1 citizen science project (Lintott et al. 2008), al¬ 
most every galaxy in the SDSS spectroscopic catalogue has 
a visual morphology classification. From the votes of the 
‘citizen scientists’, debiased from redshift-dependent resolu¬ 
tion effects that may blur galactic features (such as spiral 
arms) in SDSS images, galaxies are classified as ellipticals 
or spirals (clockwise, anti-clockwise and edge on) if either 
category receives at least 80% of the votes. Objects that do 
not reach the 80% threshold are designed as uncertain, in¬ 
cluding a fraction of merging systems. According to Darg 
et al. (2010), a visual merger identification can be consid¬ 
ered reliable when it receives more than 40% of the votes. 
Based on these criteria, we label every object in our sam¬ 
ple as early-type, late-type, uncertain (excluding mergers) or 
merger. 


2.3 Observational biases 

Although our physical characterization is intended to be as 
objective as possible, it is not completely free from observa¬ 
tional biases, some of them inherent to the selection function 
of the spectroscopic SDSS. 

For instance, the minimum fibre separation in the detec¬ 
tor implies that galaxies in the spectroscopic catalogue can 
only be as close ~ 55 arcsec in the sky (physical distances of 
~ 22.5 and ~ 74 kpc at our redshift limits of 2 = 0.02 and 
0.07, respectively). This might introduce some bias in the 
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Figure 3. Comparison between model SSFR (Brinchmann et al. 
2004) and our observational proxies. Left panels correspond to the 
SSFR within the fibre, and right panels show the total SSFR. Top 
and bottom panels show Ho equivalent width and {u — r) colour, 
respectively. Contours in blue, red, green, and orange correspond 
to the area that encloses 90% of the star-forming, passive, in¬ 
termediate, and agn subsamples. Points correspond to a random 
selection of 1500 target galaxies, and they are similarly coloured 
according to the subsample they belong to. 

distance to the fifth neighbour, and it poses a serious prob¬ 
lem for the identification of close pairs (see e.g. Behroozi 
et al. 2015, and references therein). 

In addition, most of the quantities we consider are based 
on astrometric and/or photometric measurements, but aper¬ 
ture effects may be an issue for spectroscopic observables 
such as the equivalent width of the Ha line or the 03N2 
ratio. The aperture of the SDSS fibre (3 arcsec in diameter) 
corresponds to 1 — 4 kpc at z = 0.02 — 0.07, and therefore 
these measurements may not always be representative of the 
whole galaxy. 

Aperture bias has been thoroughly discussed in the lit¬ 
erature (e.g. Brinchmann et al. 2004; Iglesias-Paramo et al. 
2013), and models that take into account the photometry 
outside the fibre have been proposed in order to correct the 
spectroscopic SSFR from aperture effects. In Figure 3 we 
compare the direct measurements of EW(Hq') and (m — r) 
with the aperture-corrected SSFR inferred from the Brinch¬ 
mann et al. (2004) prescription, as implemented in the MPA- 
JHU pipeline^. 

Equivalent width is tightly correlated with the model 
SSFR in the area covered by the fibre (top left panel), at 
least as far as star-forming and intermediate galaxies are 
concerned. For passive galaxies, the Ha absorption line is ac¬ 
tually more prominent in systems with higher SSFR due to a 
larger contribution of A stars to their optical spectrum. The 
aperture correction proposed by Brinchmann et al. (2004) 
does not change these general trends, but it broadens the 
correlation with the extrapolated SSFR, as shown on the 

® http://www.mpa-garching.mpg.de/SDSS/DR7/ 


top right panel, {u — r) colour, on the other hand, is an 
excellent proxy for the model SSFR over the whole galaxy 
(bottom right panel), and it is well correlated with the the¬ 
oretical estimate inside the fibre, even for passive galaxies 
(bottom left panel). 

The distribution of star-forming, intermediate, and pas¬ 
sive galaxies in Figure 3 clearly shows that our sample se¬ 
lection criteria are roughly equivalent to fixed thresholds in 
SSFR, colour, or equivalent width. The latter arises from 
the minimum signal-to-noise imposed, which is dominated 
by the contrast (i.e. the equivalent width rather than the 
absolute luminosity) of the H/3 line with respect to the 
underlying continuum. Our S/N > 2 criterion, and thus 
our definition of the star-forming subsample, is in practi¬ 
cal terms very similar to EW(Ha)>10 A, (u — r) < 2.3, or 
SSFR> 3 X 10“^^ which may not be statistically rep¬ 

resentative of the overall galaxy population and introduce a 
significant source of bias. 

Passive galaxies typically populate the range (u — r) > 
2.3 and SSFR < yr~^. The intermediate subsam¬ 

ple is presumably dominated by emission-line objects with 
EW(Ha) < 10 A, most likely weakly star forming galaxies, 
but, given the low signal-to-noise ratio, it also includes an 
uncertain fraction of absorption-line systems and AGN. We 
have made no attempt to correct for such effects. 

Finally, we have also assumed that the value of 03N2 
within the spectroscopic fibre is representative of the whole 
galaxy. This measurement does not account for radial abun¬ 
dance gradients, and (assuming the fibre is placed at the 
centre of the galaxy) it would be biased high for large spi¬ 
rals. The precise impact of aperture bias on our results con¬ 
cerning gas-phase metallicity is hard to predict, but the ob¬ 
served values of the radial abundance gradients in CALIFA 
disk galaxies suggest a variation of the order of ~ —0.2 dex 
between 0.3 and 2 effective radii (Sanchez et al. 2014). 


3 RESULTS 

In order to study the physical mechanisms that regulate 
star formation, discriminating between ‘nature’ and ‘nur¬ 
ture’ processes, we present in this section a thorough charac¬ 
terization of the distribution and cross-correlations between 
stellar mass, environment, close companions, morphology, 
gas-phase metallicity, and specific star formation rate, using 
the observational proxies discussed in Section 2.2. 

3.1 Star-forming galaxies 

Our main results for the subsample of star-forming galax¬ 
ies are summarized in Figure 4, where each column shows 
histograms of Mr, Rsfr, Ri/R^, galaxy Zoo morphological 
classification, 03N2, equivalent width of the Ha line, and 
[u — r) colour. In each row, different lines represent adap¬ 
tive bins in the corresponding physical property, containing 
2000 galaxies each. Numbers within each panel indicate the 
Spearman’s rank correlation coefficient between both vari¬ 
ables. 

When galaxies are binned by mass (top row), there are 
no apparent trends in the distributions of R 5 /f and RifR^. 
We find (as others before, e.g. Gomez et al. 2003) that star¬ 
forming galaxies tend to live in the field, and very few such 
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Darker color implies lower value bin 



Figure 4. From left to right, different columns show the distribution of mass (Mr), environment (Rs/f), interactions {R\/R^), mor¬ 
phology (based on the Galaxy Zoo visual classification), metallicity (03N2), and specific star formation rate, traced by EW(Ho) and 
[u — r) colour, for our subsample of star-forming galaxies. Different lines correspond to adaptive bins in another of these properties (one 
per row) containing 2000 galaxies each. Darker colours correspond to lower value bins, i.e. imply higher mass, denser environments, low 
values of R^fR^, higher metallicity, lower EW, and bluer colour, respectively. 





objects are found in the densest environments. Furthermore, 
we show here that this statement is not only valid in a 
qualitative sense; the distribution oi R^/f is independent on 
galaxy mass as long as the galaxy Is classified as star-forming 
according to our criteria (see figure 5 for a comparison to 
the intermediate and passive populations). 

In contrast, the well-known relation between metallicity 
and luminosity (Lequeux et al. 1979) is clearly visible in the 
distribution of 03N2, with more luminous galaxies (darkest 
lines) displaying systematically lower values of 03N2 (higher 


metallicities). There is a very weak trend in the distribution 
of equivalent widths (top row, sixth panel), in the sense of 
a more extended tail (i.e. a larger fraction of galaxies with 
higher SSFR, such as blue compact dwarfs and HII galax¬ 
ies) in low-mass systems. The trend Is clearer when consid¬ 
ering (n — r) colour (rightmost panel), although the values 
of the rank correlation coefficients indicate that there is a 
weaker correlation between luminosity and colour or equiv¬ 
alent width than between luminosity and metallicity. 

When galaxies are classified in terms of the environmen- 
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tal density (second row), we do not find any difference in the 
luminosity function of star-forming galaxies nor the distri¬ 
bution of their Hq equivalent widths. Nevertheless, we can 
observe that galaxies in the densest environments tend to 
be slightly more metal-rich (darkest lines in the second row, 
fifth column), in agreement with recent results (e.g. Hughes 
et al. 2013), and display marginally redder colours. 

Concerning the role of close companions, the distribu¬ 
tion of R 1 /R 2 (plotted on the panels in the third column) 
is roughly uniform, as expected for a random Poisson distri¬ 
bution in two dimensions. However, there is a clear excess 
at < 0.1 that we interpret as an indicator of in¬ 

teractions with the nearest neighbour. The only difference 
between these galaxies (represented by the darkest lines in 
the third row) and the rest of the star-forming population 
is the lack of such close pairs in dense environments. Since 
galaxies move at larger relative speeds in clusters than in 
the field, the typical type of two-body interaction is an im¬ 
pulsive encounter (‘harassment’) rather than the formation 
of a close pair and subsequent merger (see e.g. Boselli & 
Gavazzi 2006). There seems to be a very mild enhancement 
of the SSFR for the darkest bin (Ri/Ri < 0.077), observed 
in both Ha and (u — r) colour, but its statistical significance 
is rather weak. 

The fifth row is binned in gas-phase metallicity, traced 
by the 03N2 line ratio. One can see, as in the reciprocal pan¬ 
els in first and second rows, both the luminosity-metallicity 
relation and (to a much lesser extent) the slight deviation 
of the most metallic systems towards denser environments. 
A clear correlation exists between 03N2 and SSFR, with 
lower-metallicity objects showing bluer (u — r) colour and 
higher EW(Ha). In particular, about fifty per cent of the 
objects in the lowest metallicity bin (with 03N2 > 1.25, 
i.e. less than half solar according to Perez-Montero & Con- 
tini 2009) feature values of the equivalent width higher than 
50 A. 

In the sixth row, where galaxies are binned according to 
their equivalent width, no evident correlation is found with 
galaxy density or the presence of close companions. There 
is a marginal trend with luminosity, with the highest equiv¬ 
alent widths (larger than 65 A) appearing in systematically 
low-mass systems, as well as a significant correlation with 
03N2 (albeit with a smaller correlation coefficient than the 
luminosity-metallicity relation). Similar results are obtained 
using (u — r) colour (seventh row) as an alternative proxy 
for the SSFR. No correlation is found with our nurture indi¬ 
cators (Rs/f, Ri/R|), but there is a clear trend with galaxy 
luminosity. The correlation with metallicity is even stronger 
than that observed for EW(Ha). 

Concerning morphological types, the fourth column of 
Figure 4 presents the fraction of galaxies classified as late- 
type, early-type, uncertain (excluding mergers) or mergers, 
binned in terms of each of the other proxies. Conversely, 
the fourth row shows the distribution of each galactic prop¬ 
erty, segregated by our four morphological types (one can 
no longer build adaptive bins of 2000 galaxies). 

Our star-forming subsample consists mainly on objects 
with late (43.15%) and uncertain (56.28%) morphologies, 
while early-type (0.45%) and mergers (0.12%) are barely 
present. None of our ‘nurture’ proxies (Rs/f and R 1 IR 2 ) 
seems to correlate with morphology, but we must recall that 
this refers to the star-forming subsample only. The trends 


observed for the late-type and uncertain objects show that 
the former are, on average, more massive, more metallic, red¬ 
der, and display lower values of EW(Ha). Early-type star¬ 
forming galaxies (160 objects) tend to display slightly higher 
metallicities, redder colours, and higher equivalent widths. 
There are only 43 merging systems in the star-forming sub¬ 
sample, and therefore it is difficult to conclude whether there 
is a statistically significant increase in their average equiva¬ 
lent width. 

3.2 Separating ‘ageing’ and ‘quenching’ 

As shown in Section 2.3, the signal-to-noise threshold im¬ 
posed for our star-forming subsample is roughly equivalent 
to a selection criterion EW(Hq:)> 9 A. From Figure 4, we 
have just concluded that there is only a slight dependence of 
the SSFR on galaxy luminosity, and almost no dependence 
on environmental density, when such a threshold is imposed. 

We will now address the relation between galaxy mass, 
environment, and star formation, including also the inter¬ 
mediate and passive subsamples, thus covering all possible 
values of the Ha equivalent width (with our adopted sign 
convention, negative values indicate a stellar absorption fea¬ 
ture). Once again, we first arrange all galaxies in bins (either 
in Mr or Rs/r) containing 10000 objects. Galaxies within 
each bin are then classified as star-forming, intermediate, or 
passive, as described in Section 2.2, and the respective frac¬ 
tions, as well as the conditional probability distribution of 
the physical property not used for the binning {Rs/r or Mr), 
are plotted in Figure 5. Star-forming galaxies are shown on 
the top panels, whereas the central and bottom panels dis¬ 
play intermediate and passive systems, respectively. 

Let us start by focusing on the distribution of Rsfr, 
plotted on the leftmost column of the figure. While star¬ 
forming and intermediate galaxies tend to live in the field, 
passive galaxies can be found in all kinds of environments. In 
contrast with the other two populations, they also display a 
remarkably different behaviour at high and low luminosities. 
At the low-mass end (light-shaded lines) they are mostly lo¬ 
cated in dense environments, but the distribution of Rs/f 
quickly becomes bimodal as we consider higher luminosity 
bins. High-mass passive galaxies (showing Ha in absorption 
with S/N > 2) can be found both in clusters and in the 
field, confirming the idea that there are indeed two mecha¬ 
nisms that may reduce (not necessarily ‘quench’ on a short 
time scale) the star formation activity of a given galaxy to a 
negligible level. One of these mechanisms would be related 
to galaxy mass (and, most probably, its ‘nature’), whereas 
the other, clearly related to the environmental density (i.e. 
‘nurture’), would act upon galaxies of all masses (see e.g. 
Peng et al. 2010). 

As it has been previously reported in the literature, the 
fraction of star-forming galaxies (selected by either colour or 
Ha) at fixed stellar mass is a strong function of local density, 
and there is also a clear dependence on galaxy luminosity at 
fixed Rsfr. The fraction of galaxies classified as star-forming 
according to our criteria (top panels on the right column of 
Figure 5) is indeed a monotonically decreasing function of 
both luminosity and galactic overdensity, although, contrary 
to some previous studies (cf. Tanaka et al. 2004), we do not 
find evidence for a sharp transition at any particular value 
of either quantity. 
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Figure 5. Distribution (left panels) of iJs/f (environment, left) and Mr (mass, right) for the star-forming (blue), the intermediate{gTeen) 
and the passive{red} subsamples. The fraction of the three subsamples (right panels) is also represented against R^fr and Mr - Line shading 
correspond to bins in the other variable and are fixed to contain 10000 galaxies. Darker colours correspond to lower value bins, i.e., more 
massive objects and denser environments. Black solid line in left panels correspond to the distribution of the whole sample. 


Moreover, we would like to argue that the distinction 
between passive and intermediate galaxies is also a critical 
factor in order to discriminate between ‘ageing’ and ‘quench¬ 
ing’ mechanisms. In low-density environments, most of the 
galaxies that are not star-forming fall in the intermediate 
category, i.e. the fraction of truly passive systems (Ha un¬ 
ambiguously detected in absorption) is systematically lower 
than that of intermediate objects for all luminosities. As we 
will discuss later in more detail, the average SSFR of the 
field population is a decreasing function of stellar mass, but 
it seldom becomes negligible, especially for late-type galax¬ 
ies; it merely tends to drop more and more often below the 
selection threshold defining the star-forming class. 

As noticed by several authors (e.g. Baldry et al. 2006; 
Peng et al. 2010), the effects of mass and environment on 
the fraction of star-forming galaxies are fully separable, and 
the curves with different shading on both of the top pan¬ 
els of the right column are parallel to each other. However, 
these objects may gradually move towards the intermediate 
population due to ‘ageing’, staying there for a long time, 
or quickly migrate through the intermediate class towards 
the passive group due to a sharp quenching of their star 
formation activity. 

If only a mass-independent quenching mechanism were 
at work, the fraction of passive galaxies would increase at 
the expense of both star-forming and intermediate objects 


in equal proportion, and therefore the ratio between their 
fractions would stay constant. In contrast, one can readily 
observe that the effect of the environment on the fraction 
of intermediate galaxies is not independent on mass. The 
curves on the middle panels of the right column in Figure 5 
are not parallel to each other, and they bear relatively little 
relation to the fraction of star-forming galaxies on the top 
panels. 

Our results show that, as the local density increases, 
a certain fraction of the star-forming population becomes 
intermediate and another undetermined fraction becomes 
passive. At the same time, a fraction of the intermediate 
galaxies turns passive, but they are approximately balanced 
by a similar number of star-forming galaxies moving into 
the intermediate realm. Incidentally, the final fraction of in¬ 
termediate galaxies turns out to be rather insensitive to the 
environmental density (middle right panel). 

This process is further examined in Figure 6, where we 
show the distribution of SSFR, traced by both colour and 
equivalent width, in different bins of mass and environment. 
We confirm the finding of Balogh et al. (2004) that the mean 
(u — r) colour of the ‘blue’ population (a mixture of interme¬ 
diate and star-forming galaxies) becomes progressively red¬ 
der with increasing galactic overdensity. A similar trend can 
be seen in EW(Ha), whose distribution is plotted on the left 
panels. As the density increases, the distribution becomes 
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Figure 6. Normalized distributions of EW(Ha) (left) and (u — r) 
(right) for different luminosity bins (brightest galaxies are plotted 
on the top panels). Line shading indicates environmental density, 
with darker colours corresponding to the lowest values of Rs/r 
(densest environments). The binning in Rs/f has been computed 
from the whole target sample (10000 galaxies per bin), and it is 
the same in all panels. 

more asymmetric, with a tail extending into the intermedi¬ 
ate region, and its peak moves slightly towards lower values 
of the equivalent width. While these effects are consistent 
with a mild ‘ageing’ associated to denser environments, the 
dominant mechanism at play is a fast quenching of the star 
formation rate, responsible for quickly driving the objects 
towards the peak associated to the passive population. 

However, the role of intermediate galaxies is not entirely 
clear in this context. Are they transition objects, caught 
during the short interval of the transformation from star- 


formtng into passive after their SFR has suddenly dropped 
to zero, or will they live there for a long time, because their 
SFR has just decreased slightly (e.g. due to a shortage of 
external gas supply)? 

In order to discern ‘ageing’ from ‘quenching’ (i.e. the 
time derivative of the star formation rate), it is interesting 
to compare two tracers of the SSFR that are sensitive to 
different time scales. As recently shown by Schawlnski et al. 
(2014), early- and late-type galaxies occupy very different 
regions of the NUV—u — r colour-colour diagram, with the 
former displaying systematically redder (NUV—m) colours, 
even for similar values of (u — r). Since UV emission probes 
star formation on shorter time scales (of the order of 10^ — 
10® yr) than optical colours (10® — 10® yr), they concluded 
that early-type galaxies in the ‘green valley’ have undergone 
a rapid end to star formation, whereas late-type galaxies 
with similar (u — r) colours have experienced at most a slow 
decline in their star formation activity. 

Rather than UV colours, we will use here the equivalent 
width of the Ha line to probe the SSFR on scales of the order 
of ~ 10^ yr. Roughly speaking, the intensity of the Ha line 
is proportional to the mass of O and B stars, whereas the 
continuum traces the total stellar mass. Thus, EW(Ha)~ 
Mr/M, where Mr denotes the stellar mass created during 
the last 10 Myr. On the other hand, {u — r) is sensitive to 
the presence of A stars, and it varies on scales of the order 
of 300 Myr. To first order, one may consider that (m — r) ~ 
Ms.s/M, although this is only a coarse approximation (e.g. 
in the absence of dust extinction, no galaxy can be redder 
than a single stellar population with the current age of the 
universe). 

If the star formation rate of a given galaxy evolved 
slowly [ip/ip > 300 Myr), there should be a one-to-one 
relation between its Ha equivalent width and its (m — r) 
colour, since both would be tracing exactly the same quan¬ 
tity; SSFR= ip/M. We will refer to this relation (the location 
of all galaxies with smooth star formation histories on the 
EW-colour diagram) as the ‘ageing’ sequence. 

If star formation were suddenly shut off, nebular Ha 
emission would drop to negligible values on 10-Myr scales, 
but the galaxy would still retain some memory of its orig¬ 
inal blue colour. In fact, there would be a large amount 
of A stars and, therefore, recently quenched galaxies are 
expected to feature stronger absorption lines (and bluer 
colours) than a genuinely old stellar population (Dressier 
& Gunn 1983). The higher the SSFR before the quench¬ 
ing event, the stronger the Ha absorption and the bluer 
{u — r) after 10 Myr. During the subsequent Gyr, A stars 
would gradually die; Ha absorption will become weaker, and 
colours will become redder, until the galaxy eventually joins 
the red extreme of the ‘ageing’ sequence. Galaxies in this 
phase would occupy a completely different region of the EW- 
colour diagram, forming a well-defined line where EW(Ha) 
is a monotonically increasing function of (u — r) that we will 
term the ‘quenched’ sequence. 

The {u — r) colour - EW(Ha) diagram of our sample of 
SDSS galaxies is plotted in Figure 7, segregated by morpho¬ 
logical type and environment, and coloured according to the 
object luminosity. Although the ticks and labels in the y-axis 
indicate the actual value of EW(Ha), the plotted quantity 
is in fact log(EW/A-|-2) in order to show both emission and 
absorption lines on a common logarithmic scale. 
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Figure 7. Galactic colour-equivalent width diagram. Top row corresponds to the whole galactic sample. Intermediate and bottom rows 
correspond to galaxies in low (field) and high density environments respectively. From left to rigth (and in different columns), contours in 
black, blue, red, gray, and orange correspond to the area that encloses 90% of the total galaxy sample, early-type, late-type, uncertain, and 
merger subsamples. Points correspond to a random selection of target galaxies. We draw the same number of objects in the three rows 
coloured according to their mass (Mr) in three different bins: red, objects with Mr < —21.5; green objects with —21.5 < Mr < —19.5 ; 
and blue, objects with Mr > —19.5. 
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At the lowest overdensities, late-type galaxies seem to 
be well described by a single ‘ageing’ sequence (i.e. 'ip/ip > 
300 Myr). The vast majority of them are forming stars at 
some level, but they span a broad range of SSFRs, going 
from very active to very quiescent systems, and we do not 
detect any sharp transition whatsoever between the star- 
forming, intermediate, and passive populations. 

Field early-type galaxies appear at the end of this se¬ 
quence. Some of them display weak Ha emission, but almost 
none of them shows clear signatures (blue colours and strong 
Ha absorption) of a recent quenching event. Most of these 
systems are consistent with very low values of the SSFR 
(not necessarily of the instantaneous SFR) and a very old 
stellar population, in the sense of the mass-weighted average 
defined in equation (1). 

Galaxies with uncertain morphology are distributed 
along the whole sequence, but it is interesting to note that 
they reach the most extreme values at both ends (in fact, 
this category accounts for more than 80% of the galaxies 
with EW(Ha) > 65 A). From Figure 7, it becomes evident 
that galaxies showing the strongest quenching signatures 
also tend to display a disturbed appearance. Interestingly, 
field galaxies unambiguously classified as ongoing mergers 
can only be found at the blue extreme. 

In the densest environments, the ‘ageing’ sequence de¬ 
fined by late-type galaxies moves toward redder colours and 
lower equivalent widths. There is a general displacement of 
the sequence, consistent with the results plotted in Figure 6, 
in the sense that galaxies of a given luminosity tend to be 
slightly older in dense environments. However, the overall 
sequence stays roughly invariant (if anything, perhaps even 
narrower than in the field), with very few objects showing 
obvious post-quenching signatures. Therefore, we conclude 
that late-type galaxies in dense environments have system¬ 
atically lower SSFR and older stellar populations given their 
mass, but their star formation rate has not undergone any 
sudden change. 

In contrast, many early-type galaxies in dense environ¬ 
ments are arranged in a tight disposition with stronger Ha 
absorption in bluer systems, consistent with the quenching 
scenario. Furthermore, if we focus on galaxies in different 
luminosity bins, one can see them all over the ‘quenched’ 
sequence, tracing the time elapsed since star formation 
stopped, but fainter objects, which typically started with 
higher SSFR, reach bluer colours and deeper absorption fea¬ 
tures than more massive systems. 

The ‘ageing’ and ‘quenched’ sequences are also evident 
in galaxies with uncertain morphology, which constitute (as 
in the field) more than fifty per cent of the population. They 
seem to reach bluer colours in the quenched sequence, but 
they do not particularly stand out among the ageing popu¬ 
lation. A small number of mergers seem to be associated to 
objects with enhanced star formation, but most of them are 
located at the red extreme of the quenching sequence. 

3.3 The SSFR-luminosity relation 

There is ample observational evidence that there exits 
a tight relation between the stellar mass and instanta¬ 
neous/specific SFR of star-forming galaxies, often referred 
to as the ‘galactic main sequence’ (see e.g. Noeske et al. 2007; 
Speagle et al. 2014, and references therein). At low masses. 


the star formation rate is roughly proportional to the stellar 
mass, and therefore all star-forming galaxies display simi¬ 
lar SSFRs, almost independent on luminosity. The situation 
at high masses is still debated, but there is growing con¬ 
sensus that massive star-forming galaxies have consistently 
lower SSFR (albeit similar or higher SFR) than their smaller 
counterparts, especially in the local universe (e.g. Whitaker 
et al. 2014; Ilbert et al. 2014; Tasca et al. 2014). 

The SSFR of our galaxy sample, traced by either (u — 
r) colour or EW(Ha), is plotted in Eigure 8 as a function 
of galactic luminosity Mr, segregated by morphology and 
environment as in Eigure 7. In the EW(Ha) —M,. plane (top 
panels), galaxies are colour-coded by (u — r), whereas in the 
colour-magnitude diagram (bottom panels), red, green, and 
blue correspond to different bins in equivalent width. 

Since galaxies of early-type morphologies are known to 
display consistently redder colours than late-type spirals (see 
the colour-magnitude diagrams on the bottom panels; cf. 
Schawinski et al. 2014), it is not surprising that they also 
feature lower equivalent widths, often (yet not always) show¬ 
ing the Ha line in absorption (top panel). Nevertheless, some 
early-type galaxies, especially in the field, display obvious 
signs of star formation, both in EW(Ha) and in {u — r). As 
noted above, they seem to represent the low-SSFR end of 
the ‘ageing sequence’ in Figure 7. Figure 8 suggests that, in 
the absence of environmental effects (i.e. in the filed) there 
is a single relation between SSFR and Mr- 

We would strongly advocate for including intermedi¬ 
ate and passive galaxies, of any morphological type, in the 
‘galactic main sequence’. As noted by Whitaker et al. (2014), 
the shape of the sequence does not change by more than 
0.1 dex when the whole galaxy population is considered. 
Combining our Figures 7 and 8, we do not find observa¬ 
tional support for any quenching mechanism acting on field 
galaxies. The most massive systems tend to have formed 
most of their stars a long time ago and have arranged them 
into an early morphological type, but there seems to be a 
continuous variation over the full range of SSFR. 

In particular, our results strongly disagree with the two- 
population scenario, where some galaxies form stars with 
a given, mass-independent SSFR, until they are suddenly 
quenched (e.g. Peng et al. 2010). The apparent trend of a 
constant SSFR (as in e.g. Figure 4) is entirely due to the 
selection of star-forming systems. As shown in Figure 3, our 
subsamples of star-forming, intermediate, and passive galax¬ 
ies roughly correspond to fairly well defined cuts in colour 
or equivalent width. In fact, they are statistically represen¬ 
tative of the traditional ‘blue cloud’, ‘green valley’, and ‘red 
sequence’ regions selected from the colour-magnitude dia¬ 
gram. We argue that any such classification of the galaxy 
population introduces a serious bias on the average (S)SFR 
as a function of galaxy mass. 

The situation in dense environments is completely dif¬ 
ferent. The stellar population of late-type galaxies becomes, 
on average, slightly older, and some of these objects do in¬ 
deed belong to the passive class. At a given luminosity, the 
distribution of equivalent widths seems to broaden, reaching 
significantly lower values of EW(Ha), whereas the (u — r) 
colours tend to become somewhat redder, consistent with 
the ‘ageing’ effect discussed in the context of Figures 6 and 7. 
The fraction of elliptical galaxies increases dramatically with 
respect to the field, and many of them display signatures 
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Figure 8. Galactic SSFR-luminosity relation, using Ha equivalent width (top) and {u — r) colour (bottom) as proxies for the SSFR. 
Figure schema is the same as in Figure 7. The colours of the points correspond to: Top - three different bins in (u — r) colour: red, 
objects with {u — r) > 2.5; green objects with 1.5 < {u — r) < —2.5 ; and blue, objects with (u — r) < 1.5. Bottom - three different bins 
in EW(Ha): red, objects with EW(Ha)< 0.0; green objects with 0.0 <EW(Ha)< 10.0 ; and blue, objects with EW(Ha)> 10.0. 
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that are consistent with recent quenching. Both statements 
are particularly evident in the case of low-mass galaxies, 
which typically display the highest SSFRs at low overdensl- 
ties, and therefore the effect of the environment can be more 
clearly noticed as a sharp discontinuity between the galaxy 
properties of the ‘ageing’ and ‘quenched’ populations. 

Galaxies with merging signatures display a clearly bl- 
modal distribution. In the field, they feature the highest 
SSFR, especially as far as the EW(Ha) is concerned, whereas 
in dense environments they are systematically located at 
the red extreme, although it is not obvious whether they 
are more likely to belong to the ‘ageing’ or the ‘quenched’ 
sequence (cf. Figure 7). 


3.4 The SSFR-metallicity relation of star-forming 
galaxies 

Our results so far are consistent with a single ‘galactic age¬ 
ing sequence’ where most galaxies spend the largest part 
of their lives, including both passive and star-forming sys¬ 
tems, with early- or late-type morphologies. In addition to 
the correlations discussed in the present work, many other 
relations between the relative properties (colours, equiva¬ 
lent width, spectral shape, metallicity, gas fraction, mass- 
weighted mean age of the stellar population...) of a given 
galaxy have previously been reported (see e.g. Ascasibar & 
Sanchez Almeida 2011; Ascasibar et al. 2014), and therefore 
any of these parameters can be used to describe the cur¬ 
rent state of the galaxy in terms of its location along the 
sequence. 

Here we will attempt to use the gas-phase oxygen abun¬ 
dance, traced by the 03N2 line ratio. As shown in Fig¬ 
ure 4, the SSFR (measured in terms of colour or equivalent 
width) presents a higher correlation coefficient with 03N2 
than it does with the absolute magnitude Mr. On theoretical 
grounds, the gas-phase metallicity is expected to be closely 
related to the Integrated star formation history of the galaxy, 
not only for the closed-box model (Searle & Sargent 1972), 
but for more elaborate scenarios including inflow and out¬ 
flow of gas (e.g. Edmunds 1990), even on resolved scales 
(Ascasibar et al. 2014). The main disadvantage, of course, is 
that enough signal-to-noise ratio is necessary in order to es¬ 
timate the 03N2 line ratio, and therefore we have restrained 
the discussion to the star-forming subsample. 

The observed relation between the SSFR and 03N2 
of star-forming galaxies, illustrated in Figure 9, is indeed 
tighter than the relation between SSFR and luminosity dis¬ 
cussed in the previous section. As it was previously reported 
when discussing Figure 4, the trend is clearer when using 
{u — r) colour as a proxy for SSFR. The SSFR-metalllcity 
relation is fully compatible with ‘ageing’ dominating the 
evolution of most star-forming galaxies, displaying almost 
exclusively late-type and uncertain morphologies. Interest¬ 
ingly, we do see that the few mergers present in the sample 
do not only show, on average, higher SSFR, but also lower 
chemical abundance, than late-type or uncertain systems, 
hinting that they may be associated to the accretion of low- 
metallicity gas towards the main galaxy. In all cases, galaxies 
with high SSFR (i.e. younger stellar populations) are metal- 
poor, and vice versa. Additional work would be required 
in order to investigate whether the observed correlation be¬ 


tween SSFR and metallicity is not only stronger, but indeed 
more fundamental than the correlation with stellar mass. 

Regarding denser environments (recall here that we are 
missing passive and intermediate galaxies) we find that there 
is an evident lack of low-metallicity objects that could be 
explained either by an accelerated ageing (Induced by e.g 
a reduction of the gas accretion rate with respect to the 
field) or by unidentified quenching mechanisms (e.g. strip¬ 
ping and/or rapid consumption) completely removing the 
gas from our systems and sending them to the quenched 
sequence. In the latter case, it might be possible to find a 
quenching signature in the relation between SSFR and the 
mass-weighted stellar metallicity. 


4 DISCUSSION 

In this section, we would like to propose a scenario where a 
combination of both nature and nurture processes is respon¬ 
sible for regulating the SSFR of a given galaxy. As noted 
in the introduction, we distinguish between smooth varia¬ 
tions of the instantaneous/specific SFR (ageing), a sudden 
interruption of the star formation activity (quenching), and 
a transient episode of star formation that is sufficiently in¬ 
tense as to decrease the mass-weighted age of the underlying 
stellar population (rejuvenation burst). Let us now briefly 
explore some possible physical mechanisms that may be re¬ 
sponsible for these qualitatively different behaviours. 


4.1 Ageing 

In low-density environments, we argue that most galaxies 
should be expected to start their lives in a ‘chemically- 
primitive’ (gas-rich, metal-poor, young stellar population, 
and high SSFR) state. As they turn their gas into stars and 
enrich it with metals, they would gradually move towards a 
‘chemically-evolved’ (gas-poor, metal-rich, old stellar popu¬ 
lation, and low SSFR) state. Let us stress once again that 
this ‘ageing’ process is unavoidable unless the star formation 
rate increases exponentially in time. While possible, such 
sudden changes in the SSFR (see the subsection on ‘rejuve¬ 
nation bursts’ below) do not seem to dominate the overall 
evolution of typical field galaxies, although they may signif¬ 
icantly perturb it for a limited amount of time and/or for a 
limited fraction of objects. 

Considering the average properties of the galaxy popu¬ 
lation, unaffected by environmental factors, this work firmly 
establishes the existence of an ‘ageing sequence’; i.e. a fairly 
tight relation between the average EW(Ha) and {u — r) 
colour of unperturbed galaxies (extensive to 03N2, at least 
for the star-forming subsample). This sequence is not meant 
to represent chronological evolution, but a particular region 
of the parameter space where most galaxies are found at the 
present time (and there is, in fact, no reason to expect that 
it should be invariant with cosmic epoch). It merely reflects 
the relation that must exist between the physical properties 
of ‘young’ and ‘old’ galaxies if their star formation rate has 
varied slowly with time according to some smooth function 
ij){t), whose analytical description lies well beyond the scope 
of the present discussion. 

Quite remarkably, dwarf galaxies (and, in particular. 
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Figure 9. Galactic SSFR-metallicity relation, using Hq equivalent width (top) and {u — r) colour (bottom) as proxies for the SSFR, 
and 03N2 for the gas-phase oxygen abundance. Figure schema is the same as in Figure 7. The colours of the points correspond to: both, 
top and bottom - three different bins in mass {Mr) : red, objects with Mr < —21.5; green objects with —21.5 < Mr < —19.5 ; and blue, 
objects with Mr > —19.5. 
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HII galaxies) tend to be located closer to the primitive ex¬ 
treme of the sequence, whereas more luminous systems tend 
to be, on average, more chemically evolved. The physical 
origin of the relation between the total mass and the mass- 
weighted age of the stellar population (or the gas-phase 
metallicity) is an open question. Ageing is related to the 
speed at which the gas is accreted and turned into stars, 
which may be affected by several ‘nature’ and ‘nurture’ pro¬ 
cesses. The observed dependence on the environmental con¬ 
ditions may thus help us constrain the dominant physical 
mechanisms at work. 

Galaxies in the field are consistent with a single trend, 
without any obvious discontinuity associated to the stel¬ 
lar mass or SSFR. For ‘ageing’ galaxies at fixed luminos¬ 
ity, there is a mild shift towards more chemically-evolved 
systems (higher metallicities and lower SSFR) in dense en¬ 
vironments. Since both the stellar mass and the galaxy over¬ 
density are related to the mass of the host dark matter halo, 
this shift might be associated with ‘nature’ processes driven 
by the relative efhciencies of gas infall, cooling, and con¬ 
sumption as a function of halo mass. In addition, ‘nurture’ 
processes, such as the strangulation of cold gas supply, may 
play a role in the accelerated ‘ageing’ in dense environments. 

As a final remark, let us note that galaxies with a 
smooth star formation history may be subject to morpho¬ 
logical evolution, triggered by both internal and external 
mechanisms, and they may also be temporarily offset from 
the ‘ageing sequence’ due to discrete quenching or rejuve¬ 
nation events. As long as they remain in the field, they will 
continue accreting gas and forming stars. After some time, 
they will return to the ‘ageing sequence’ and recover their 
original morphology, which may be most likely late-type or 
uncertain. 

The time scales of such recovery will depend on the 
mass of the object, the availability of gas, and intensity of 
the event under consideration. We argue that field elliptical 
galaxies, located at the end of the ‘ageing sequence’, are 
representative of the last stages of secular galactic evolution. 
We do not observe any convincing evidence that they have 
recently interrupted their star formation activity due to a 
merger event. Rather, we propose that these galaxies, being 
more massive, already had a lower SSFR previous to such 
event, and therefore they kept their post-merger elliptical 
morphology up to the present day. In this scenario, early- 
type morphologies would be a consequence, rather than the 
cause, of a reduced SSFR. 

4.2 Quenching 

One of the main results of the present work is that we only 
find unambiguous evidence for quenching in dense environ¬ 
ments. In these cases, star formation does not evolve grad¬ 
ually until the gas reservoir is exhausted. On a very short 
time scale, these galaxies are removed from the ‘ageing se¬ 
quence’ by shutting down star formation to a level that Ha 
is clearly detected in absorption. Then, they evolve over the 
lifetime of A stars along a ‘quenched sequence’ in the colour- 
equivalent width diagram, until both branches merge at the 
‘red-and-dead’ extreme. 

The critical role played by the galactic overdensity hints 
that quenching is directly related to the mass of the dark 
matter halo. It seems reasonable to expect that, for star for¬ 


mation to stop, it is necessary to prevent further accretion 
of cold gas. This would be the reason why most galaxies in 
the field (in particular, low-mass galaxies) are almost invari¬ 
ably found among the ‘ageing’ population. Even if a discrete 
event might momentarily quench their star formation activ¬ 
ity, it would be quickly resumed as a consequence of efficient 
gas accretion and cooling. 

Nevertheless, the very existence of a ‘quenched se¬ 
quence’ implies that the physical process(es) responsible for 
environmental quenching must be very fast. Preventing the 
accretion of new cold gas is a necessary but not sufficient 
condition; the galaxies’ reservoir must be almost entirely 
consumed on a time scale much shorter than the lifetime 
of A stars. Other ‘nurture’ processes, such as ram-pressure 
compression and stripping, or tidal interactions and mergers 
with other galaxies, could contribute to quench star forma¬ 
tion in dense environments, probably leaving a morphologi¬ 
cal signature. From the top panel of Figure 8, we do indeed 
find that most of the low-mass galaxies featuring Ha absorp¬ 
tion lines have early-type or uncertain morphologies. Some 
of them are classified as spirals, though, and it is difficult 
to assess whether the events (most likely nurture-related) 
responsible for the morphological transformation have also 
played a dominant role in the quenching of star formation, 
or, by the contrary, the observed morphology is a conse¬ 
quence of quenching, analogously to the scenario we pro¬ 
posed for field elliptical galaxies. 

4.3 Rejuvenation bursts 

Both ‘ageing’ and ‘quenching’ drive our systems towards 
a more ‘chemically evolved’ state. In particular cases (ar¬ 
guably much more common in the early universe), it is pos¬ 
sible that some galaxies evolve towards more ‘chemically 
primitive’ states for a limited period of time. For instance, 
accretion of pristine gas well above the star formation rate 
would reduce the stellar-to-gas fraction and dilute the chem¬ 
ical composition. It is also likely that such event drives a 
burst of star formation that would increase the metal abun¬ 
dance and stellar fraction, as well as the SSFR, yielding a 
higher equivalent width for the duration of the burst as well 
as higher luminosities and bluer colours for a longer time. In 
extreme cases (regardless of gas accretion), a star-formation 
burst would decrease the mass-weighted average age of the 
stellar population. This is our dehnition of ‘rejuvenation’. 

Assuming that these processes leave an imprint on 
galaxy morphology, one can assess their statistical rele¬ 
vance on galaxy evolution by studying the correlations be¬ 
tween colour, equivalent width, luminosity, and metallicity 
for galaxies displaying disturbed morphologies. Our results 
clearly show that visually-identihed mergers display system¬ 
atically higher EW(Ha) and lower (u — r) for a given Mr 
than galaxies with late-type morphology (or even those with 
uncertain morphology) for low-mass galaxies in low-density 
environments, i.e. those where a significant gas reservoir is 
expected. We interpret this as evidence for merger-induced 
rejuvenation. 

However, most galaxies with high SSFR have uncer¬ 
tain morphologies, which could well be associated to a re¬ 
cent merger, but they could also be due to intrinsic ‘nature’ 
processes (e.g. stochastic star formation), especially in low- 
mass galaxies. We can only ascertain that galaxies with ir- 


© 2015 RAS, MNRAS 000, 1-18 


16 J. Casado et al. 


regular morphology have, on average, slightly higher SSFR 
than late-types, although the effect seems to be milder than 
the enhancement observed in ongoing mergers (which rep¬ 
resent a very small fraction of the galactic population). We 
argue that ‘rejuvenation’ episodes contribute to the scatter 
observed in both, the ‘ageing sequence’ and the relations be¬ 
tween SSFR, metallicity, and luminosity. Nevertheless, these 
processes are not the dominant mechanism regulating nei¬ 
ther the star formation nor the chemical evolution in galax¬ 
ies. 

It is possible, though, that they are important in partic¬ 
ular types of objects, such as e.g. HII galaxies. The presence 
of an underlying host with an older stellar population (see 
e.g. Papaderos et al. 1996; Amorfn et al. 2009; Janowiecki 
& Salzer 2014) hints that these objects are good candidates 
for a ‘rejuvenation’ episode. However, it is not entirely clear 
whether the mass-weighted age of the stellar population, de¬ 
fined in equation (1), has decreased as a consequence of a 
recent star formation burst or the SSFR has been a smooth 
function of time throughout most of their history. Recent 
observations do indeed suggest that low-mass star-forming 
galaxies, including blue compact dwarfs, have formed 90 per 
cent of their stellar mass between 0.5 and 1.8 Gyr ago, con¬ 
sistent with the expectation from the SSFR-luminosity (or 
SSFS-metallicity) relation of ‘ageing’ galaxies (Rodrfguez- 
Munoz et al. 2015). 

Since EW(Ha) is sensitive to star formation ^ 10-Myr 
scales, we argue that close pairs or merger signatures should 
be clearly visible in the majority systems in the interaction- 
induced burst scenario. Therefore, our results support the 
idea that HII galaxies are simply the ‘chemically primitive’ 
extreme of the ‘ageing sequence’, but we do not find evidence 
that their SSFR has increased drastically on very short time 
scales. 


5 CONCLUSIONS 

We have carried out an analysis of the mechanisms driv¬ 
ing and regulating star formation in the local universe. To 
do so we have selected a sample of galaxies from the SDSS 
DR7 spectroscopic catalogue, applying restrictions in appar¬ 
ent magnitude, redshift, and position in the sky, and then 
classified them as star-forming (35425 objects), agn (8099 
objects), passive (19266 objects), or intermediate (19826 ob¬ 
jects) according to the observed intensity and signal-to-noise 
ratio of some of their optical spectral lines. 

For the star-forming subsample, we have studied the 
distribution of absolute magnitudes in the r-band {Mr), the 
normalized distance to the fifth nearest neighbour (Rs/f), 
the ratio of squared projected distances of the first two 
neighbours {R^/R^), the morphologies derived from the 
galaxy Zoo visual classification project, the 03N2 metal¬ 
licity indicator, the equivalent width of the Ha line, and 
the {u — r) colour. Although we find evidence that mergers 
may temporarily increase the SSFR, we conclude that most 
(non-merging) star-forming systems evolve through internal 
processes, i.e. ‘nature’, as we do not find, on average, any 
statistically-significant signature of the interaction-induced 
scenario. We do find, however, that star-forming galaxies 
tend to be found in the field, and we show that the distri¬ 
bution of environmental density is in fact independent on 


other galaxy properties, including luminosity. As far as we 
are aware, this is the first time such result has been reported 
in the literature. 

Our analysis also shows that the restriction in signal- 
to-noise used to define our star-forming subsample, roughly 
equivalent to a fixed threshold in EW(Ha) or (it —r), creates 
an artificial bimodality when comparing star-forming and 
passive galaxies. Including the intermediate population, we 
arrive at the following conclusions: 

(i) In low-density environments, most present-day galax¬ 
ies are distributed along a relatively narrow ‘ageing se¬ 
quence’ in the EW(Ha)-(it — r) plane. At the ‘chemically 
young’ extreme, galaxies have high SSFR (blue colours and 
high equivalent widths), low metallicity (high 03N2), and 
late-type or uncertain morphology. We argue that, at the 
‘chemically old’ extreme, galaxies tend to display elliptical 
morphologies because the SSFR is insufficient to rebuild a 
significant disk. Brighter galaxies tend to be more chemi¬ 
cally evolved, resulting in clear correlations between SSFR 
and luminosity or metallicity, but all of them are consistent 
with a smoothly-varying SSFR driven by secular (‘nature’) 
processes, with little evidence for recent ‘quenching’ or ‘re¬ 
juvenation’ episodes. 

(ii) In dense environments, we detect a ‘quenched se¬ 
quence’ in the EW(Ha)-(it —r) plane, consistent with a very 
rapid truncation of the star formation activity. Most of these 
objects display early-type or uncertain morphologies, but it 
is difficult to establish a causal connection between galaxy 
interactions and the drop in SSFR. 

We thus propose a scenario where ‘nature’ is more im¬ 
portant than ‘nurture’ in regulating star formation in galax¬ 
ies. Starting their lives as ‘chemically primitive’ objects, they 
gradually turn their gas into stars with a slowly-evolving 
SSFR and evolve towards ‘chemically old’ systems, lead¬ 
ing to the observed relation between EW(Ha) and (it — r) 
colour. Apart from this unavoidable ‘ageing’ process, galax¬ 
ies are also susceptible to ‘quenching’ (in dense environ¬ 
ments) and/or ‘rejuvenation’ episodes (associated to ex¬ 
treme bursts of star formation). The ability to return to 
the ‘ageing sequence’ after such events is linked to the ex¬ 
ternal gas supply. Galaxies will become red and dead either 
when they totally exhaust their gas reservoir and stop form¬ 
ing stars (end of the ‘ageing sequence’) or when they suf¬ 
fer a quenching event in a dense environment and they are 
no longer capable of accreting gas (‘quenched sequence’). 
Mergers and galaxy-galaxy interactions may temporarily af¬ 
fect the instantaneous SSFR, but they merely seem to add 
statistical fluctuations to the main relation. 
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